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Table 1: List of stellar characteristics of the stars sample: coordinates J2000, equivalent 
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intrinsic unaccuracies of the TMGS coordinates and the extreme crowding in the field there 
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ABSTRACT 

This paper is the second part of Garzon et al. (1997) in which we presented 
an outline of the analysis of 60 spectra from a follow-up program to the Two 
Micron Galactic Survey (TMGS, Garzon et al. 1993) project in the £ = 27°, 
b = 0° area. In this second part, we present a more detailed explanation of the 
analysis as well a library of the spectra for more complete information for each 
of the 60 stars, and further discussions on the implications for the structure of 
the Galaxy. 

This region contains a prominent excess in the flux distribution and star 
counts previously observed in several spectral ranges, notably in the TMGS. 
More than 50% of the spectra of the stars detected with rrix < 5.0 mag, within 
a very high confidence level, correspond to stars of luminosity class I (Garzon 
et al. 1997), and a significant proportion of the remainder are very late giants 
which must also be rapidly evolving. We make the case, using all the available 
evidence, that we are observing a region at the nearer end of the Galactic bar, 
where the Scutum spiral arm breaks away, and that this is powerful evidence for 
the presence of the bar. Regions of this type can form due to the concentrations 
of shocked gas where a galactic bar meets a spiral arm, as is observed at the 
ends of the bars of many face-on external galaxies. Alternative explanations do 
not give nearly such a satisfactory account of the observations. 

Equivalent spectroscopical analysis should also be performed at £ — —22°, 
the candidate position for the other tip of the bar. The space localization of one 
and, a fortiori, of both ends of the bar allows us to infer its orientation. If the 
second region is also confirmed to be a powerful star formation region this would 
imply a position angle for the bar of ~ 75° with respect to the Sun-Galactic 
centre line. This geometry is indeed compatible with the range of distances 
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that we have obtained for the star-forming region at i = 27° from spectroscopic 
parallaxes. However, the angle is different from that given by other authors for 
the bar and this, we think, is because they refer to the triaxial bulge and not to 
the bar as detected here. 

Subject headings: stars: formation - Galaxy: stellar contents - Galaxy: structure 
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1. Introduction 

This paper is the second part of Garzon et al. (1997) (hereafter G97) in which we 
presented an outline of the analysis of 60 spectra from a follow-up program to the TMGS 
project in the i = 27°, b = 0° area. This analysis was based on the luminosity classification 
of stars using the Call triplet (Jones et al. 1984; Diaz et al. 1989), and a discussion of this 
region was presented which strengthens the identification of this zone as a star formation 
region at the receding tip of the stellar bar. 

In this second paper, we present the spectra, offer a more detailed explanation of 
the analysis, and supply a more complete information for each of the 60 stars, extended 
bibliography, and further discussions of the physical morphology. 

2. Observations and data reduction 

As explained in G97, 70 spectra of stars situated in the region around £ = 27°, 6 = 0° 
selected from the TMGS with m K < 5.0 mag were taken with the Isaac Newton Telescope 
in La Palma (Canary Islands, Spain). From these, we had to reject ten which could not 
be identified as visible counterparts of TMGS sources. After cross-correlating the original 
TMGS source positions with their Guide Star Catalogue counterparts and then calculating 
the errors in position, we inferred that TMGS catalog has a positional error box of around 4 
arcsec in right ascension by 7 arcsec in declination. The larger error in declination compared 
with that in right ascension is due to the form of detector array used in the survey (Garzon 
et al. 1993). This is added to the pointing error of the INT telescope itself and implies 
that there is a non-negligible risk of taking the spectrum of the wrong star in such a very 
crowded area of the sky. By taking spectra of all the bright stars within the error box, we 
could decide post hoc which candidate was the TMGS object with m K < 5.0 mag from the 
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spectral type and features and /-magnitude estimate. In the ten cases mentioned above, we 
could not find any good candidate so these were rejected. 

In Table [5] and Fig. |2], coordinates of all the stars in our sample from the TMGS 
database are given. Investigators who may wish to use them should be aware that, due to 
the error of the TMGS coordinates described above, and the extremely crowded nature of 
the field there is often more than one visible candidate for the infrared source. As we did 
not use more precise coordinates for a given source, we have not needed to sharpen the 
precision of these positions. 

The extraction of the spectra, the sky subtraction and the wavelength calibration 
were performed using the IRAF tasks apextract and onedspec after performing the 
bias subtraction and the flat fielding with ccdred. The flux was normalized rather than 
calibrated. 

3. Luminosity classification 

Spectra in the zone of the IR Ca II triplet at 8498.02 A, 8542.09 A, and 8662.14 A, 
present in stars of spectral types later than F5, were examined. The triplet was chosen as 
an optimum differentiator between stars of different luminosity classes. In earlier types, the 
Paschen hydrogen lines P13 (8665 A), P15 (8545 A) and P16 (8502 A) severely contaminate 
the spectral region of interest. For spectral types later than M3 or M4, TiO absorption 
bands mask the Ca II triplet, making its features more difficult to measure, as we will see 
later. 

Measurement of the equivalent width of the Ca II triplet lines was carried out according 
to the prescriptions of Diaz et al. (1989) when TiO-band contamination is not present. 
A continuum /™ nt is defined by a linear fit to the median value of two chosen side-bands 
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centered at 8455 A and 8850 A, respectively, with a width of 15 A and the equivalent width 
was obtained by integrating: 

where I\ is the measured intensity. 

We also followed Diaz et al. (1989) in defining 



EW = J¥ ° + W (2) 

8542 A 8667 A v ' 

as the indicator of luminosity class. The shortest-wavelength line of the triplet was not 
used because it is too feeble compared with the neighboring features and yields larger errors 
in the precise measurement of the total Call triplet. Selected reference stars were used to 
check our EWs and these were the same, within the error margin, as those given by Diaz et 
al. (1989) for the same objects. 

The Call triplet has been used for many years as an indicator of luminosity class (e.g., 
Merrill 1934). More recently Jones et al. (1984), and Diaz et al. (1989) have calibrated 
empirically the relationship between the equivalent width of the Ca II triplet and the 
luminosity class for spectral types ranging from F5 to M3. These authors found some 
dependence of EW on both metallicity and temperature, but these effects were much weaker 
than the dependence on surface gravity, especially for supergiants. We have followed their 
results and adopted their criteria in assigning luminosity classes from the measured EWs. 
From Diaz et al. (1989), if metallicity dependence is assumed negligible, then 

log g = 7.75- 0.65EW. (3) 

On this basis, a source is taken to be a supergiant (log g less than approximately 1.75) 
if EW > 9 A, independently of metallicity and temperature. Idiart et al. (1997) have 
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demonstrated the presence of correlations of EW with metallicity and temperature, but 
only for old stars (> 1 Gyr), and therefore not for supergiants. There are also other possible 
luminosity-class indicators. In the wavelength range we used, the existence of the CN band 
(7916-7941 A) is a characteristic of supergiants (MacConnell et al. 1992), but the Ca II 
triplet is less affected by noise. 

Thirty-eight stars in the sample were not contaminated by the TiO band, and we used 
the method described above to obtain their EWs in this subsample (see Table H). 

The remaining 22 stars are of later spectral types, and the presence of TiO bands 
partially masks the triplet lines (see Fig. || c), d)). For these objects we have developed 
an empirical method which permits the measurement of EW where the Ca II lines are not 
completely masked by the TiO band. This method uses the depth of the lines instead of 
the EW. 

The depth of the lines is calculated according to 



mm 



(J V ;A- 15 A < A' < A + 15A) 



^pth A = 1 — s (4) 

where J^ ont2 is the value of a new continuum interpolated between two points that are not 
affected by the TiO absorption: 8432 A and 8844 A Q. 

One potential problem is that the convolution of the profiles of the spectral lines with 
the instrumental profile would affect the depth of the line but this effect is in practice 
negligible compared with the other uncertantities in the method. The convolution with the 
instrumental profile broadens the line and decreases the depth when the number of pixels 



x In fact, a bandwidth of 15 A between 8417 A and 8432 A is selected for which we 
calculate the median value, and the same is done for the second value of the continuum 
between 8844 A and 8859 A. 
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defining the line is small. In the lines without TiO contamination, the pixels defining a line 
cover a sufficient wavelength range for the effect of the convolution to be negligible. In the 
lines with TiO contamination, the depth is the sum of the TiO band depth and the jump 
between the TiO band and the peak of the line. If the unblended portion of the line covers 
a restricted wavelength range, this jump will be affected by the convolution but, as it is 
small compared with the TiO depth, the effect of the convolution is again small. The effect 
of convolution would indeed be serious if we were to use Gaussian fits on masked portions 
of the lines instead of the depth as the measure of equivalent width. 

To get the EW as a function of the depth of the deepest line, 8542 A, it is necessary to 
use uncontaminated lines. Expression (B) is used for the first 38 stars, those with unmasked 
Call lines, and for these a relationship is then calibrated between depth 854 2 and the EW 
calculated using eq. (Q). The relationship obtained, employing a linear fit, was: 

EW = -0.1 + 17.8 x depth 8542 , (5) 

which is plotted, together with the data, in Figure 00. The use of this expression to evaluate 
the luminosity class using the criterion of Diaz et al. (1989) is a fair approximation. The 
fit in equation (||) is a good approximation but in the latest spectral types other effects 
may intervene. The equivalent width of the Ca II triplet is not strongly dependent on 
spectral type, although this dependence is a little more noticeable for M-type stars since 

The fit to a second-degree polynomial was EW = -7.1 +45.0 xdepth 8542 - 25.8 xdepth^ 
which is very close to the linear fit and yields the same values, within the error limits. We 
emphasize that in G97 we used a second-degree fit, but with a small error, which gave rise 
to EW values too low by 0.2 or 0.3 A. This led to our classifying two stars as giants instead 
of supergiants. This improved result only strengthens our previous conclusions about the 
proportion of supergiants in the sample. 
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the ionization equlibrium in these objects shifts somewhat from Ca II to Ca I (Cohen 1978). 
One consequence is that this line-depth technique may slightly underestimate the number 
of supergiants. That is, the number of supergiants among stars later than M3 may well be 
in fact larger than that inferred using the EW > 9 criterion applied to equation (|5]). Since 
our purpose is to calculate the fraction of supergiants to determine whether the region 
studied is a star-formation region, this underestimate could change our conclusion only if 
the answer were negative, but this is not the will be shown below. 

Even employing the line-depth technique, we had to reject two stars of the sample 
for which the TiO band gave an unacceptable blend, so we finished with 58 stars — 38 
earlier than M3, and 20 later than M3 — with acceptable EWs. The two stars rejected show 
spectra with no clearly unblended Call features; the Ca il+TiO blends show a maximum 
depth ° from which we may estimate a maximum EW. In both cases, this maximum 

^ 8542 A 

value of EW is ^10.2 A, i.e. EW < 10.2, so it could not be decided whether they are giants 
(EW < 9.0) or supergiants (EW > 9.0). In Table ||, we do show the data for these two stars 
with EW < 10.2 (stars 29 and 37) but these are not included in either the statistics of the 
EW or in the distance calculation of paper G97. 

The final results of the luminosity classification were shown in Fig. 1 of G97, in the 
form of a histogram of EW frequencies. Those EW values used for the figure are shown in 
Table |6| for each star. 

The ratio of SGs — those with EW > 9 A — to the total number of observed objects is 
most striking. The number of SGs is in fact 38 out of 58 (66%). In a binomial distribution 
the root mean square of the distribution is a = J~n x p x (1 — p), where n is the total 
number of stars. Here, n = 58 stars, p = 0.66, so a = 3.6 stars. The proportion of 
supergiants is then 66 ± 6% (1 a). An error should be added to allow for the possible 
mistaken classification of giants and supergiants, and also to take into account the pointing 
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error of the telescope in a crowded field which gives a slight chance of taking the spectrum 
of a wrong star. However, if the "wrong" star turned out to be a late star with sufficient 
luminosity to give a spectrum of bright star, this would still be indication of late type 
giants and supergiants in the zone. These sources of error are small except in the M-type 
classification, in which the number of supergiants should be even greater that given by our 
criterion, as discussed before. These considerations lead to the conclusion that the fraction 
of supergiants is well over 50% with a very high confidence level. 

4. Spectral classification 

The spectra of the 60 stars belonging to both subsamples are presented in Fig. ||, after 
removing the intrinsic slopes and those due to reddening, in order to show the qualitative 
differences in the spectra between different kinds of stars. The spectral type classification is 
carried out for every star by qualitatively comparing the features of our spectra with those 
of standard stars. The spectral classification is shown in Table 

The coolest stars — M3 and later on — were compared with those in Bessel (1991), 
Schulte-Ladbeck (1988), and Barbieri et al. (1981), by inspecting the depth of the TiO 
band absorption. In Fig |6] c) and d), it can be seen that there is a growing depression 
at 8432 A and 8844 A from M3 to M9. Since the features of these sources permit us to 
differentiate spectral types with an interval of 0.1 or 0.2 spectral types, the classification is 
quite accurate for these stars. 

The earlier-type stars were compared with those in Torres-Dodgen & Weaver (1993), 
and Schulte-Ladbeck (1988), mainly by inspecting lines such as Mg I 8807 A, Fe I 8612 A, 
Fe I— Ti I 8468 A, O I-Fe I 8446 A (some of these features are also dependent on luminosity 
class or metallicity) and Paschen hydrogen lines for those earlier than G2. In this case the 
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enors in the classification could be as large as half a spectral type, except perhaps for those 
containing Paschen lines. 

A histogram of the frequency of spectral types in our sample is shown in Figs. 2 and 3 
of G97. As expected, most of the stars are very red and with intrinsic luminosity (see also 
Calbet et al. 1996a) since TMGS is more sensitive to redder stars. The predominant type 
among the supergiants is K. Supergiants of very late types are known to exist in very small 
numbers; we have found a significant number of new examples here. This result is subject 
to the bias that our method of predicting the EWs gives where the TiO band affects the 
Ca II lines; it tends to underestimate the EWs, thereby reducing the apparent fraction of 
SGs in the coolest (M) class. 

5. Stellar bar 

The conclusion arised from this analysis is that the observed region contains intense 
star formation of recent origin, since supergiants and bright giants are necesarily young 
stars. There is a high proportion of young stars, and their high spatial concentration shows 
that the star formation has taken place in the zone observed. Any associated early B-type 
or O-type stars, typical of this kind of cluster, are not observed in this case since their 
emission in K is too low to reach m K < 5 mag for objects far from the solar neighborhood, 
and the emission in the visible is very strongly affected by dust extinction. 

This star formation region does not belong to the disc nor the spiral arms (G97); 
nor it is likely belonging to the ring, as exposed in G97 and Hammersley et al. (1994) 
(hereafter H94). A further convincing reason for excluding a ring, even an elliptical one, 
besides the ones in those papers, is that it would be prominent in other TMGS regions 
closer to the center than £=27°, and that the luminosity function of the stellar excess 
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would be invariant, neither of which would be in agreement with the observations. The 
TMGS star counts after subtraction of a model disc (Wainscoat et al. 1992) and bulge 
(Lopez-Corredoira et al. 1997) are zero in off-plane directions, showing that those disk and 
bulge models are good fits to the observations; there is, however, an excess of stars in the 
plane (\b\ < 2°). Figure [| shows a sharp decrease in the star counts in the center, and also 
that the luminosity function of the stellar excess at £=27° is different from that in the other 
regions (the ratio between stars with tuk < 9.0 mag and mx < 5.0 mag is very different). 
These considerations lead to the conclusion that a ring-type structure, even elongated and 
eccentric, gives a poor fit to the observations. 

A hole in the extinction was suggested as a possible cause of the excess at £=27° by 
Jones et al. (1981) but this idea was criticized by Kawara et al. (1982) who observed 
an invariant (H-K) colour for the stars across this region. Also, the TMGS star counts 
at £=27° in the plane show an excess of over 100 stars per square degree up to m K = 5 
and this cannot be satisfactorily explained by assuming that the extinction is zero in that 
region, according to our calculations. Furthermore, CCD VRI photometry of the region 
shows clearly that the line-of-sight extinction is greater than the standard value of 0.62 mag 
kpc -1 in V (Mahoney 1999, Hammersley et al. 1998). 

The most likely explanation is the presence of a giant star formation region which is 
associated with the receding tip of the stellar Galactic bar (see Fig. |3[). If we take the 
region to be in the center of the star formation region at the end of a barQ, it can be 
used not only as a qualitative demonstration of the bar's existence but also as a means to 

3 Its center is at £ ~ 27°, b ~ 0°. More precisely, Viallefond et al. (1980) places the region 
with center at £ = 27.5°. The region is extended over several degrees but this is the point 
where the star counts show a maximum which we believe correspond to the center. 
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estimate its orientation. The remaining stars in Fig. [5] can be explained as bar stars, with 
a prominent peak at both ends due to the interaction with spiral arms. Such regions can 
form due to the concentrations of shocked gas where the Galactic stellar bar meets a spiral 
arm, as is observed at the ends of the bars of face-on external galaxies (Sandage & Bedke 
1994). An excess of extincion at negative latitudes (Calbet el al. 19966) is also explained 
by the stellar-bar hypothesis. 

This is not the first time that a Galactic bar has been claimed to be discovered but our 
arguments in favour of it are new and different from those of other authors. De Vaucouleurs 
(1964, 1970) first suggested that the Galaxy might be barred in an attempt to explain the 
observed non-circular gas orbits. Since then, many types of observational evidence have 
been accumulated that support this hypothesis (see Blitz 1993; Blitz et al. 1993; Kuijken 
1996a; Gerhard et al. 1997). The axial asymmetry of the inner Galaxy is detected by star 
counts (Nakada et al. 1991; Weinberg 1992; Whitelock et al. 1992; Stanek et al. 1994, 1996; 
Wozniak & Stanek 1996; Nikolaev & Weinberg 1997), or by surface photometry at different 
wavelengths (Blitz & Spergel 1991; Weiland et al. 1994; Dwek et al. 1995, Sevenster 1996), 
stellar population studies in Baade's window (Ng et al. 1996), micro- lensing (Paczyhski et 
al. 1994) and studies of the extinction using lensed stars (Stanek 1995), analysis of internal 
motions of the gas (Peters 1975; Liszt & Burton 1980; Yuan 1984; Gerhard 1996). Models 
including a bar (e.g., Binney et al. 1991; Zhao 1996) can explain these observational 
features. 

Not all of these observations necessarily imply the existence of a bar; some could be 
explained by a triaxial structure in the inner Galaxy. Whether this structure is a triaxial 
bulge, or whether a bar or both features are present is still a matter of controversy (Kuijken 
19966; Ng 1997), and there are reasons for believing that many authors who refer to a bar 
are in fact referring to the triaxial bulge. 
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In our previous analysis of the Two Micron Galactic Survey (TMGS) database 
(Garzon et al. 1993) we described evidence in favor of the presence of a triaxial bulge 
( Lopez- Corredoira et al. 1997) with radial scale length ~2.2 kpc, making an angle 12° 
with respect to the Sun-Galactic center line in the first quadrant, in strong agreement with 
the angle given by many of other authors (Binney et al. 1991; Weinberg 1992; Dwek et 
al. 1995; Binney et al. 1997; Stanek et al. 1997; Nikolaev & Weinberg 1997) whose bar 
angle (around 20 degrees) is quite close to this. The bar of this paper must be a different 
structure from the triaxial bulge since the inclination of the bulge is not sufficient for a 
bar tip to reach I = 27° nor for a dust lane to reach £ = —19° at the other side of the 
bar, as found by Calbet et al. (19966). The triaxial bulge reported by Lopez- Corredoira et 
al. (1997) does not yield sufficient star counts in the plane at positive Galactic latitudes, 
taking into account the disk and the extinction, to explain the observations. 

There is another infrared peak at negative Galactic longitudes: at £ = —22°, which is 
very similar to the one analyzed herefj} If we assume that the other end of the bar is at 
£ = —22°, the orientation of the bar can be derived. In fact, it is a simple trigonometrical 
problem when we take the bar to be rectilinear and with equal length from each end to the 
center. 

Applying the sine rule (see Fig. [|), 

Lp -Ro 

sin 27° ~ sin(180° - 27° - a) ' U 

4 See COBE-DIKBE data (Boggess et al. 1992) or near-infrared catalogs which cover the 
Galactic plane at positive and negative Galactic longitudes. See also H94, and Calbet et al. 
(19966). 
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Lp _ Rp 
sin 22° ~ sin(a - 22°) ' 

Hence, with Ro = 7.86 kpc (Lopez-Corredoira et al. 1997), we deduce that a 
and the length of the bar from tip to tip is 2L = 7.4 kpc. 

Therefore, we suspect that the authors which are giving a much lower angle are in fact 
analysing the angle of the triaxial bulge, or an average angle of the composition of both 
bulge and bar, instead the bar angle. The triaxial bulge cannot be causing the features seen 
at £ = 27°. 

We can also derive geometrically the distance to the tip at £ = 27° as being 7.8 kpc. 
This distance is compatible with that derived in G97 (in table |6] are the distances for each 
star), taking into account all the uncertainties in the calculation. 

6. Conclusions 

This paper, together with its first part (G97) makes a spectroscopic analysis of the 
brightest stars in an infrared selected sample of objects close to the Galactic plane at 
I = 27° showing a strikingly high fraction of supergiants, characteristic of a strong star 
formation region. 

We argue, using all the available evidence, that this region is situated at one end of 
the Galactic bar, where the Scutum spiral arm breaks away, and that the presence of the 
region is itself a powerful evidence for the presence of the bar. None of the alternative 
possibilities (arm, disk, bulge or ring components) is capable of explaining the observations 
in a satisfactory manner. 

The detected presence of a similar concentration of near-IR sources in the plane at 



(7) 
= 75.6° 
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£ = —22° should, on this hypothesis, indicate the other end of the bar. To confirm this 
requires a similar spectroscopic campaign to that reported here. A rectilinear bar between 
these two points makes an angle of 75° with the Sun-Galactic centre line, and has a total 
length of 7.4 kpc. The distance this implies to the star-forming zone analyzed here is 
consistent with the estimated spectroscopic parallaxes of the stars whose spectra we have 
analyzed in the first part of the paper (G97), giving a self-consistent picture of the bar. 

It is important to distinguish this bar from the triaxial bulge of the Galaxy (Lopez- 
Corredoira et al. 1997). Many of the previous claims to describe a Galactic bar are more 
likely to refer to the triaxial bulge and not to the bar. 

The Isaac Newton Telescope is operated on the island of La Palma by the Royal 
Greenwich Observatory in the Spanish Observatorio del Roque de Los Muchachos of 
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PB94-1107, and PB97-0219 of the Spanish DGICYT. 
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FIGURE CAPTIONS 

Figure 1: Equatorial coordinates (J2000) of the stars whose spectra are analyzed in this 
paper. Small circles stand for giants and big circles stand for supergiants. 

Figure 2: Comparison between EW measured in unmasked Call lines using (|2[) and the 
subsample having the EW obtained via fl5|). 

Figure 3: TMGS star counts as functions of Galactic longitudes in the Galactic plane 
after subtraction of the Wainscoat et al. (1992) model disc and the Lopez- Corredoira 
et al. (1997) model bulge. 

Figure 4: A bar can explain a number of observations: the excess of extinction at negative 
Galactic latitudes due to a dust lane (Calbet et al. 1996b), and the star formation 
regions at each end of the bar. 

Figure 5: Geometry of the bar sustaining an angle a with respect to the Sun (S)-Galactic 
center (C) line. 

Figure 6: Plots of the spectra from our sample in the range 7900A- 8900A. The intrinsic 
slopes and those due to reddening have been removed. Numbers for the stars in the 
right hand side correspond to those given in the first column of the Table || and the 
spectral type and luminosity class of the star. 



